Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN WF$t style file v2.2) 



Foreground contributions to the Cosmic Microwave 
Background 



On 



T. Wibig u and A. W. Wolfendale 2 

1 Experimental Physics Dept., University of Lodz; The Andrzej Soltan Institute For Nuclear Studies, 
Cosmic Ray Lab., Lodz, Uniwersytecka 5, POB 447, Lodz 1; Poland 

■ 2 Physics Department, University of Durham, Durham DHI 3LE, UK. 

o ■ 

£h 1 2 February 2008 
t-H ' ABSTRACT 

<n : 

A detailed search has been made for evidence of re sidual foreground contr ibutions 
<N ■ to the Cosmic Microwave Background (CMB) WMAP llBennett et alJl2003alblh a map 

(nominally) cleaned for foreground already l)Tegmark et alJ 120031) . We find positive 
results in that various features relate to Galactic properties. For example, on the largest 
| angular scales we find significant differences between the power in the fluctuations for 

. positive and negative Galactic latitudes and between the four Galactic Quadrants. 

' There are also differences between the power spectrum at latitudes within 10° of the 

t^J- \ Plane and at higher latitudes. The mean temperature shows similar variations. An 

explanation in terms of Galactic effects seems inescapable. 

In an effort to find the origin of these Galactic-style effects we have examined 
the evidence from Galactic gamma rays, specifically from the EGR ET instrument 
iHunter et al.lll997|h We are mindful that the CMB maps examined iTegmark et alJ 
O ' l2003[) have already been 'cleaned' (for cosmic ray and other effects) in a rather complex 

+3 . way but, in our view, the cleaning has left potentially serious 'contaminations'. A 

correlation is found between gamma ray intensities and the CMB and other cosmic 
\ ray indicators. For example, regions of the Galaxy having (line of sight) steep cosmic 

. ray energy spectra have low mean CMB temperatures and the important Loop I edge 

J^j ' region, where the cosmic ray intensity is high, has a high mean temperature. 

rS | Most of the large scale Galactic asymmetries (eg North, South difference and 

■ Quadrant variations) have analogues in cosmic ray asymmetries and also in some 

other Galactic properties, such as the column density of gas. Thus, it is possible 
to hypothesize on direct cosmic ray-induced contributions, although it may be that 
cosmic rays are simply the indicators of Galactic 'conditions' which are influencing 
the residual CMB fluctuations. 

Irrespective of the actual cause of the correlations we have endeavoured to extrap- 
olate to the situation where the residual foreground is minimised. The effect on the 
usually derived cosmological properties is briefly examined. The least that can be said 
is that the 'error' in some of these properties has been underestimated. 
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1 INTRODUCTION 

It is well known that although the Galactic foreground 
to the CMB measurements should be small, it is not 
zero. Some 15 years ago w e feandav fc Wolfendale! |l990; 
iBandav fc Wolfendale 1 1991 drew attention to the fact that, 
despite the contribution to the 'CMB' spectrum from elec- 
tron bremsstrahlung falling rapidly with frequency, the con- 
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tribution from dust rises rapidly and there is a frequency 
region - from 30 - 100 GHz - where their sum has nearly 
the same 2.7K spectrum. The derived value of AT/T was 
about 2.5 x 10~ 6 , ie AT is about 7/iK. T he angular scale was 
8° relating to the Tenerife experiment dWatson et allll98gl) 
(beam width 8.8° FWHM and the beam throw 8.2°), and 
the value related to a 'quiet' region of the Galaxy, 'quiet' in 
the sense of having low radio and dust signals. The corre- 
sponding value for b — 90° was close to this value. It was 
claimed, further, that the minimum foreground (synchrotron 
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plus dust) would be at a frequency of 70 GHz for the 'DC 
signal', ie the absolute intensity as distinct from the fluctu- 
ations in the foreground temperature. An important param- 
eter is the degree of flatness of the foreground temperature 
vs frequency. Taking a straight line between 20 and 90 GHz 
as datum, the maximum d ownward deviation was a factor 
3.0. Actual measurements fee Oliveira-Costa. et"al"ll2004h 
now enable the DC level to be determined. WMAP data 
at the five frequencies (22.8 to 93.5 GHz) do indeed show a 
minimum signal at about 60 GHz and the form of the depen- 
dence of temperature on frequency is even flatter, the factor 
being only 2.0. The absolute temperature level is somewhat 
higher, even, than 'our' early prediction. 

Taking our predicted temperature fluctuations, as dis- 
tinct from the DC level, at AT ~ 7/j.K for the minimum, the 
average over much of the Galaxy, but away from the Galac- 
tic Plane, would be of ord er 20uK. Inspection o f the power 
spectrum for the WMAP feennett et alJl2003lJ) shows that 
for 8° (£ ~ 36) the average power is ~ 1000 /j,K 2 . 'Our' 400 
fj,K 2 is therefore a significant fraction. There is thus a prima 
facie case for a serious contribution. 

Correction of the WMAP for th e Galactic foreground 
has been made by the WMAP group feennett et alJl2003 J) 
who have carried out a procedure to generate, not a pub- 
lished map but a power spectrum. Many workers have used 
this power spectrum to derive important c osmological pa- 
rame ters. Here we will use a cleaned map llTegmark et alJ 
2003) the cleaning of which was carried out by a more so- 
phisticated technique, specifically the individual frequency 
sets o f data were weighted n ot just according to frequency 
(as in lBennett et all <l2003cf0 but also by harmonic {£) in 
the power spectrum analysis. 

An important point to be mentioned at this stage is the 
fact that the well-known cosmic ray (electron synchrotron) 
contribution to the foreground at low frequencies is cor- 
related with that of dust at the highest frequencies. This 
is a consequence of the well-known strong correlation be- 
tween the synchrotron and far infra-red radiation emission 
of galaxies, commented on by many workers ove r the years, 
including ourselves eg llChi fc WolfendaleJll990l) . Although 
there can be no certainty about the mechanism, a possibil- 
ity relates to CR-heating of the very cold dust clouds in the 
Galac tic Halo. In this resp ect we are mindful of the com- 
ment feurton fe Deullll987h that "... between a few K and 
about 30 K, the (neutral) gas is almost entirely in the form 
of molecular hydrogen, residing in dense clumps...". Inter- 
estingly, confining attention to cold clumps, their contribu- 
tion could well fall off with Galactic latitude more slowly 
than coscc \b\, because many of the low |6| clouds would 
be warmed by Galactic Plane radiation and not recorded in 
the CMB map. In fact, the ^-dependence may not be at all 
simple. 

Before continuing it is instructive to examine the var- 
ious energy densities. The CMB radiation has an energy 
density of 0.24 eVcm -3 and is universal; that in cosmic rays 
is ~ 0.5 eVcm -3 , and, as is well known, this is the same as 
that in interstellar magnetic fields, gas motion and starlight 
- all in the local region of the Galaxy. The fluctuations in 
the CMB correspond to about (1-10) xlO" 6 eV cm -3 - on 
various angular scales - and a fraction of the CR energy 
density (or of one of the other components of similar energy 
density) of order 10 -6 going into pseudo-CMB signals would 



not appear impossible. In view of the importance of the CR 
'background' in other searches, eg for cosmic neutrinos, dark 
matter particles, magnetic monopoles, it would perhaps be 
surprising if CR effects were absent. 

How much of the CR foreground has in fact been taken 
out by the 'cleaning' is open to question but we presume 
that most of the contribution has been removed and that 
the residual will vary slowly with Galactic latitude. Inso- 
far as the synchrotron radiation has a pronounced asymme- 
try about the Galact ic Plane, with a North/South (N/S) 
ratio less than unity feennett et al.ll2003d : iBurton fc Deull 
Il987fl - and all t he WMAP frequenci es show this effect 
before correction feennett et al J 12003 J) we will search for 
such an asymmetry in the cleaned map; such an asymme- 
try is also present in other Galactic phenomena, as will be 
demonstrated later. Similarly, in view of all Galactic phe- 
nomena being dependent on the Galactic Quadrant con- 
cerned, Quadrant variations will also be searched for. 

Our work is not the first to claim asym metries in the 
cleane d WMAP. Although the WMAP group feennett et alJ 
claimed that there is no evidence for non-Gaussian 
primordial fluctuations in the WMAP, and they quote many 
other workers who have come to t he same conclusion, other 
recent studies (Hans en et alJ 2004l. have shown that 'strong 
non- Gaussian features emerge when considering the North- 
ern and Southern Galactic hemispheres separate ly'. Others 
still dEricksen. et al J 120043 : lEriksen et alJl2004bl) had, pre- 
viously, found asymmetries in the distribution of the CMB 
fluctuations in the two hemispheres. Thus, the present work 
is another contribution to the debate on asymmetries in gen- 
eral, but with an independent approach. 

The cosmic ray aspect is, so far as we know, truly orig- 
inal. 

Some remarks are necessary first about the adopted 
map used throughout the analysis. This is the 'high resolu- 
tion foreground- clea ned CMB map from WMAP' given by 
iTegmark et al.H2003l) . Although the WMAP group itself has 
analysed the basic data from the WMAP instrument in great 
detail, and deri ved best-fit results for the various cosmologi- 
cal parameters ('B ennett et alJ 2003a b), it has not given the 
necessary cleaned-map data. Furthermore, its treatment of 
the five WMAP frequency maps to allow for foreground ap- 
pears to us to be rudimentary in that i'-independent statisti- 
cal weights were applied bef ore the data were combined. In 
the map adopted by us, ( Tegmark et alJl2003T) . ^-dependent 
weights were used and, although we contend that significant 
foreground remains, the data should be better. 

We estimate that the resolution of the adopted data is 
about 12' and that, although the minimum in the power 
spectrum at I ~ 350 is enhanced somewhat, the effect on 
the 'one-degree-peak' should be small. The authors state 
that the Galactic residuals seen along the Galactic Plane 
'are present in only a tiny fraction of the total sky area and 
therefore contribute little to the total power spectrum'. 



2 LARGE SCALE ASYMMETRIES 

2.1 The Parameters Under Study 

The CMB sky is characterized by the temperatures of indi- 
vidual pixels and these can be studied in a variety of ways, 
as follows: 
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i) The search for discrete excesses and deficits. For sizes 
of the order of a few pixels most will relate to extragalactic 
systems such as active galactic nuclei. We have identified 
several AGN and the VIRGO, PEGASUS 1 and PISCES 
clusters. In general, we find six times as many excesses as 
we find from searches in 'artificial universes', taken from 
the conventional CMB power spectrum. Turning to deficits 
(eg Sunyaev-Zeldovich reductions in galaxy clusters) we find 
twice as many in the genuine map as in artificial universes. 
The significance of the above is two fold; our technique of 
analysis receives a measure of support; the integrated mag- 
nitude of their contributions is small - and only affects the 
results at the highest of £ values. 

ii) An analysis of the power spectrum formed from the ba- 
sic data. The standard technique is usually adopted to give 
W{1) = £(£+l)C t /2ir, where £ is the multipole and C e is the 
angular power spectrum average derived by averaging the 
sum of the squares of the terms in the spherical harmonic. 
The expression gives the power per logarith mic intervals of 
£ and i nspection of the form of W{£), eg jBennett et alJ 
2003bK lTegmark et alj2 003). shows that much of the power 
is at large £, specifically in the region of the 'one-degree 
peak'. We expect the Galactic foreground effects to be on 
the scales of rather small £ (perhaps £ = 10) as distinct 
from the £ — 200 related to the 'one-degree peak' but all 
scales are, in principle, possible. Here, we usually study the 
^-dependence of parameters but, where undue complexity 
would arise, we give results for just £ = 10 and £ = 100. 

iii) Averaging the sky over large regions gives the mean 
CMB temperature, (AT)cmb- By definition, over the whole 
sky this will be zero, but for restricted regions the value 
will be finite, and its study can be useful. Inspection of the 
ATcmb distribution for individual pixels shows a median ex- 
cess of 80 \iK (of those which are possible) with 10% having 
AT > 200 [iK. The 80 p,K level in the standard WMAP 
maps is the 'green level', with many regions having angular 
size of order several tens of degrees. 

From what has been stated it is apparent that it will be 
valuable to look at features in the map from both the power 
- and the mean temperature - points of view. 

2.2 North-South Asymmetry 

For all Galactic radiations there is a major symmetry, with 
respect to latitude, together with a N-S asymmetry (details 
of which are given later). The former usually predominates 
but the effect of the map cleaning would be expected to 
remove most of this; the result is that the symmetry and 
asymmetry may be of not too dissimilar magnitude. A near 
cosec \b\ dependence (but bearing in mind previous remarks) 
might be expected; in fact, any disk-like emission (or ab- 
sorption) would be expected to give such a |6|-dependence, 
assuming near constant emissivity per unit volume over 
the Galactic disc. It should be remarked th at the original 
WMA P results for the individual frequencies feennett et alJ 
l2003cl) . show a pronounced cosec \b\ dependence over the 
range for which data are given: cosec|6| = 1 to 4, the slope 
of which is systematically higher in the South than the North 
(there are also subtleties in the results to which we will refer 
later). Thus, there are both symmetries and asymmetries. 
A basic question is "to what extent has the cosec \b\ de- 



pendence, and with it the largest scale Galactic foreground, 
been taken out in the 'cleaning-process' ?" 

The fact that the cosec |6| dependence, and N-S asym- 
metry in the basic WMAP data, are independent of fre- 
quency is an important result in that it suggests the fact 
mentioned in Sect. that there may be a spatial connection 
between the 'low' frequency emission (22 GHz) , which is re- 
lated to cosmic ray electrons, and the 'high' frequency emis- 
sion (93 GHz), which is affected by Galactic dust. This was 
the logic behind adding the foregroun d AT values from these 
two mechanisms in our earlier wor k. dBandav fc Wolfendalei 
ll990UBandav fc Wolfendale.lll99ll) . We start by comparing 
the power spectra for the N & S for the same latitude bands. 
We are mindful of the 'cosmic variance' which pervades the 
interpretation of all the WMAP data, but particularly at 
low values of £. Figure shows the situation. The smooth 
curves through the points for the North have been given in 
the graph for the South, for illustration purposed only. Inter- 
pretation needs to allow for the non- independence of nearby 
^-values for a particular set of data; we estimate that this 
dependence exists over a range in £ of A£ a factor 2. 

Furthermore, there is the diminishing effect on the 
South/North ratio at small \b\ and low £; for example, at 
£ = 10 (~ 18°) the ratio should fall at |6| < 10° where the 
innermost fa-values are separated by less than 20°. It does; 
from 3.1 at |6| = 10° - 20° to 1.3 for \b\ < 2°. 

The most important plot is that of the South/North ra- 
tio vs £. Figure lc shows the result summed over all latitudes; 
also shown is the mean spread for cosmic variance alone. It 
is interesting to note that the S-excess persists over virtually 
the whole range of £. Formally, it is most significant at the 
largest ^-values, a surprising result. The statistical signifi- 
cance of the excess is not easy to estimate but, allowing for 
the dependence of nearby values (within A£ ~ a factor of 
2), it is at the 1% level. 

Turning to the 6-dependence of the S-N ratio, the sig- 
nificance is far less than might be expected from the large 
excesses seen in Figure 01, because of the small solid an- 
gles for the individual b-ranges. Furthermore, the excesses 
are greatest at low £ where lack of independence of adja- 
cent bins, for small b, are most serious, as remarked already. 
However, the consequent errors are not large; inspection of 
Figures ^ and^l shows that the weighted (by solid angle) 
ratios from Figure [Ua, are 3.14 (£=5) and 2.1(^=10), to be 
compared with 2.7 and 1.7 respectively. The 'error' is thus 
only about 20%. 

Restricting attention to £ < 10 and \b\ > 10°, the de- 
pendence is less serious and something useful can be said 
Allowing for the reduced solid angles, the excesses are in the 
region of the 'upper limit' to the variance for 10° < | b\ < 60° , 
a region having 70% of the solid angle. The corresponding 
significance of the southern excess, after allowing for the fact 
that the measured ratio rises at a faster rate than the upper 
limit expected by chance, gives a formal chance probability 
of ~3%, this value relating to the low ^-values (much of the 
weight for the previous 1% came from higher I- values). 

The results on the N/S ratio suggest that there is a fi- 
nite Galactic foreground, affecting both low £ (all b) - and, 
indeed, perhaps all £. In the next section we continue to 
study the evidence for foreground effects by examining the 
data Quadrant by Quadrant (Quadrant 1 covers the range 
£ = 0° to 90°; Quadrant 2, the range 90° - 180°, and so on. 
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Quadrants 1 and 4 comprise the Inner Galaxy and Quad- 
rants 2 and 3 the Outer Galaxy). 

2.3 Quadrant-dependent effects 

There are well-known differences in all Galactic properties, 
Quadrant to Quadrant. Near the Galactic Plane, at least, 
the Inner Galaxy has higher intensities of every radiation 
than the Outer Galaxy. We start with the power spectrum. 

Figure [5^, shows the mean power vs £ for each Quad- 
rant for latitudes below 10° and Figure |!J> shows results for 
latitudes above 10°. There is little doubt that the two \b\ 
regions differ. Figures^: and[5]i show ratios with respect to 
Quadrant 1. 

It is evident that a case can be made for there being big 
variations in power, Quadrant to Quadrant, particularly at 
small £ and for |fo| < 10° these persist to high I- values. 

Considering, first, |6| < 10°, and £ = 10, the variance ex- 
pected has an upper limit ~ 4.6, so that Q3/Q1 and Q4/Q1 
are on this limit. Q2/Q1 is below it. The overall chance prob- 
ability for i < 10 is ~10%. At higher £, where the variance 
is smaller, there is an interesting and significant difference 
between the Inner Galaxy (Ql & Q4) and the Outer (Q2 
and Q3), the power in the Inner is decidedly higher. Inner 
versus Outer is different - for £ > 50 - at about the 1% level. 

Turning to |fc| > 10°, the solid angle is bigger (83% 
of the whole) but the excess is smaller. Interestingly, the 
significance of 'the effect' is similar, at least for Ql, Q2 and 
Q3 with respect to Ql, in the region of £ ~ 5. At higher £, 
the significance drops, although above £ ~ 100, the ratios 
for Q3 and Q4 are formally significantly above unity. 

In terms of the absolute magnitudes, Figures |2£ and|2f 
show comparisons between the power as a function of £ for 
the two Inner and Oute r Galaxy Quadrants a nd that for the 
standard CMB model (iBennett et al.ll2003bl) . divided by 4 
(ie each for one quarter of the solid angle). 

Some remarks are necessary about these two Figures. 
Firstly, considering the small 'wiggles', these are not signif- 
icant. Secondly, the difference between Inner Galaxy (ex- 
emplified by Ql) and Outer Galaxy (Q3) is seen. Thirdly, 
the sum of the powers for each Quadrant is nearer conven- 
tional expectation for Ql than Q3, a result due to the larger 
contribution from |6| < 10°. 

Further analysis comes from an examination of the over- 
all mean CMB temperature (as distinct from the power, 
AT 2 ) as a function of latitude. Figure [3] shows the results. 
Inspection of the figure shows a number of features 

i) There are marked differences Quadrant to Quadrant in 
mean temperature. Quantification of a lower limit to the 
significance of the variations can be made by assuming that 
the variations within one Quadrant from one latitude band 
to another are representative of the errors and deriving the 
mean displacement from zero. Considering |fe| > 10°, where 
successive latitude bands are sufficiently independent, we 
find, taking the 8 'regions' (4 Quadrants, N and S); 4 cases 
between and 1 standard deviation, 3 between 1 and 2 and 
one as high as 5 standard deviations. There is thus evidence 
for significant excursions from zero. 

In the Galactic Plane itself (|6| < 2°) there is a clear ex- 
cess temperature. The average is 27 \iK and, if it is valid to 
determine the error from the dispersion of the four sets of 



data this discrepancy from zero is at the 4 standard devia- 
tion level. Using the mean error for the other 8 sets of data 
at |b| > 10° (and five values for each set) the significance 
rises to nearly 6 standard deviations. The excess is far big- 
ger than o ne would have expecte d from the remarks about 
cleaning in iTeemark et al.l J2003D . 

Equally important are the symmetries, as follows: 

ii) Figures and [3Jj show clear symmetries in the b de- 
pendence for the two Quadrants in the Inner Galaxy (Ql 
and Q4) and those in the Outer Galaxy (Q2 and Q3). Turn- 
ing to data for the N and S separately Figures OJc) and (d), 
in Quadrants 1 and 2 there is fair symmetry between N and 
S. In Quadrants 3 and 4 there is some symmetry (but less 
than for Quadrants 1 and 2) for \b\ > 10° (it is relevant to 
point out that the giant 'hole' in the CMB, at £ ~ 300° to 
360°, b ~ 0° to —30°, is responsible for the very low mean 
temperature in Quadrant 4 for b: —2° to —20°. The sym- 
metry can be quantified by comparing the separation of the 
two values, at each |6|, with ^/2 times the modules of the 
distribution of the individual values about zero. The result 
for the ratio of the median values, (equivalent to a number 
of standard deviations) is, for Ql, Q2, Q3 and Q4; 2.1, 2.1, 
1.4 and 1.1. The symmetry is seen to be reasonably signifi- 
cant for Ql and Q2. Taken overall, N-S symmetry is present 
at the 2.4 sigma level. 

iii) There is no systematic excess temperature in the South 
as there is in the power spectrum at low I (Figure at 
least averaged over all b (but see later). 

iv) In the Outer Galaxy (Quadrant 2 and 3) there is some 
evidence for a minimum in the range |6| = 30° — 40°, ie a 
symmetrical feature. 

v) In the Inner Galaxy (Quadrants 1 and 4) there is a dip 
in the range 2° — 20° in all four cases, again a symmetrical 
feature. 

vi) In every case the mean temperature for |6| < 2° is 
higher than that for |&| = 2° — 6° and higher than that 
averaged over the |6| = 2° — 20° range. 

vii) There is a fall in (AT)cmb with increasing |6| over the 
range |6| = 6° — 45°, after which there is an increase. 

viii) Figures[3{e) and|3f) show good symmetry between the 
two Inner Quadrants (1 and 4), and between the two Outer 
Quadrants (2 and 3). 

The conclusions to be drawn at this stage from the results 
on (AT) are: 

i) There is a clear excess at \b\ < 2° averaged over all 
Quadrants and thus there is certainly some foreground left 
in this range. 

ii) Whilst there are clear Galactic features these are not of 
simple cosec b form over the whole range of |6|; nevertheless 
there is a reduction of temperature from \b\ =6° to 45°, the 
region with most of the solid angle of the sky. 

iii) There are interesting symmetries about the Galactic 
Plane. 

Comparison with the other results on 'Gaussianity', etc will 
be made later. 

The results of the mean temperature are not necessarily 
inconsistent with those for the power spectrum (Figure |2J 
because the mean temperature relates to a form of integral 
over all ^-values, whereas the power spectrum refers to the 
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square of the temperature and is, of course, ^-dependent. 
This aspect can best be studied by examining plots similar 
to Figure for the power (per unit solid angle) for different 
^-values. We remember that the power spectrum represents 
the power per unit logarithmic interval in £. 

Figure 0] shows the results for two values of I, I = 10 
and t = 100. The main feature here of relevance to symmetry 
are: 

i) Again, there is an excess nearest the Galactic Plane - in 
this case for |6| < 10°. The significance has been examined as 
for mean temperature, deriving the mean value for the other 
latitudes, and its error, and finding the difference between 
the value for |&| < 10° and this mean. There are 8 sets of data 
and the result is that for £ = 10 the average difference is 2.2 
standard deviations - ie over 5 standard deviations overall. 
For £=100, the significance is even bigger, amounting to 10 
standard deviations. For I = 100, where the interdependence 
of different |6|-bands is least, there is a strong difference 
between Inner and Outer Galaxy. In the range |6| < 10°, the 
excess in the Inner Galaxy is 2.6 times that in the Outer. 
The significance is 3.2 standard deviations. For \b\ > 10° the 
excess changes to a deficit. The ratio is 0.72 ± 0.05, ie 5.6 
standard deviations from unity. There is clearly no evidence 
for a constancy of power. 

ii) In every one of the 16 cases (for both I- values) the power 
in the |fe| = 10° — 20° is lower than would be expected on 
the basis of the adjacent points; thus a Galactic symmetry. 

iii) For £=10, in three out of four Quadrants there is a slow 
fall of power with increasing latitude beyond |6| = 20°. A 
plot of y/ (AT 2 ) vs cosec |6| in this range yields a slope of 

0.3 + 0. 04/-0. 15 (^A" 2 sr- 1 ) 1/2 for \ b \ > 20°. The mean S/N 
ratio in this region can be derived from Figure 01. It is 1.3 
(the square root of the derived value, which related to AT 2 ), 
to be compared with y/1.7 = 1.3, for all longitudes at I = 1, 
(Figure 0). An important result is that the mean slope for 
the South is greater than that for the North, although this 
result comes mainly from Quadrant 4. None of the features 
described above would be expected for 'cosmic variance'. 

iv) At £ = 100, the effects are too small, in comparison 
with the errors, for a pattern to be discerned with respect 
to |&| dependence and N-S ratio. 

The conclusions here are that, again, there is residual 
foreground near the Galactic Plane, the range |b| = 10° —20° 
appears to have been overcorrected and that, at low £, there 
is an all-latitude foreground. Furthermore, there is evidence 
for a cosec \b\ dependence of the root mean square power (at 
the 2<t level). 

Turning to the N-S asymmetry, we note that the S/N 
ratio for power, is, for £=10, and increasing \b\: 2.11, 2.49, 
1.79, 1.85, 2.25 and 0.95. leaving out the (anomalous) Quad- 
rant 4, the S/N ratios are: 2.39, 2.69, 1.53, 2.83, 1.50 and 
0.72. 

The S/N ratio therefore dips towards unity at \b\ ~ 
—60° in the second case. Inspection of the 408 MHz re- 
sults feennett et alJl2003ch . shows that there is a reduc- 
tion in ratio starting at \b\ ~ 40° reaching a minimum at 
1 6 ~ 60°. The individual frequency bands for WMAP (un- 
cleaned) show a ratio of ~1.37(+0. 13-0.04) at |6|=30° (and 
little variation at lower 6), 1.21 (+0.04 - 0.02) at |6|= 40° 
and 1.06 (+0.04-0.02)at |6|=60°, the limits encompassing 
the whole range of frequencies. It is interesting to note that 



there is no systematic change to the S/N ratio with fre- 
quency; but both in the 408 MHz and the individual WMAP 
frequencies there is a reduction in ratio at mid-latitudes; 40° 
at 408 MHz and ~40° in the individual WMAP frequency 
AT-values, to be compared with 45° - 60° in the £= 10 
WMAP power. Once again, then, there is evidence of insuf- 
ficient cleaning. 

An important aspect is a comparison of the two Quad- 
rants in the outer Galaxy, Q2 and Q3 in both (T) and power; 
we would expect greater activity in Q2 in comparison with 
Q3 because of the greater number of SNR, HII regions, etc 
there. 

Inspection of Figure[3]shows that the temperature in Q2 
is higher than that in Q3 for every latitude bin, except for 
the (anomalous) region |6| :~ 10° =20° and the previous bin. 
The excess (Q2-Q3) in AT is successively, for |6j = 0° - 2°, 
2°-6°, , > 60°: +17, +12, -7, -5, 7, 3, 3 and 19 (J.K. An excess 
is thus well founded. Insofar as much of AT comes from high 
£- values (and thus small angular sizes) , there should be little 
cross-talk between the two Quadrants. 

Surprisingly, perhaps, the same feature (Q2>Q3) does 
not pertain for the power per steradian (Figure [IJ. For £ = 
10, only for the highest latitudes (jbj > 60°) is Q2>Q3 - 
and here we would expect the nearest approach to equality. 
For |6| > 10° the ratio for Q2/Q3 is 0.56, for \b\ = 10° - 20° 
it is 0.43 and the average for |6| > 20° is 0.69. At £ = 
100, a similar situation pertains: for |6| > 10°, Q2/Q3 = 
0.93, for |6| > 10° - 20°, Q2/Q3 = 0.97 and for |6| > 20° 
Q2/Q3 = 0.94. Again, in all latitude bins, Q2/Q3M with 
the exception of \b\ > 60° where Q2/Q3 = 1. 

Thus, there is a clear difference in the power per stera- 
dian for Q2 and Q3 and one that is systematic over all lati- 
tudes (see Figure At first sight it is of the opposite sense 
to that expected but the answer may lie in the fact that the 
radial gradient of the CR intensity is ( surprisingly) s maller 
for Q3 in comparison with that for Q2 dErlvkin et alj fl996 % ). 

Turning to the S/N ratios for Q2 and Q3, overall (over 
all b) these are: Q2, 1.63 (1.04); Q3, 1.79 (1.09). The first 
value related to £ = 10 and the second to £ = 100. 

The corresponding values for b > 20° are: 1.69 (0.96) 
and 1.72 (1.02). The ratio is thus rather stable with respect 
to the limit on \b\. A different situation pertains for (T), 
where the differences (S-N) are, for Q2, +33 {iK (+16) and 
for Q3 -6 (J.K (+25 nK) the first value is for \b\ > 2° and 
that in brackets is for |6| > 20°. The fact that a S-excess is 
present for \b\ > 20° in both power and (T) is reassuring, as 
is the similar magnitude for both Quadrants. 

2.4 Comparison with the original WMAP and 
408 MHz data 

A further aspect of the 408 MHz data is the presence of 'fine 
structure' in the cosec |&|-plots' and also in the separate plots 
of AT vs cosec |6|. In every one of the six plots dWatson et all 
1988) there are consistent features, in order of diminishing 
magnitude: 

- Minima near cosec \b\ = 2.0 (ie |6| = 30°) 

- Minima at cosec \b\ = 3.15 (|6| = 18.5°) 

- Minima at cosec |6| = 1.15 {\b\ = 60.4°), all for the 
North. 

For the South there are: 



6 T. Wibig and A. W. Wolfendale 



- Minima at cosec |6| = 1.15 (|6| = 60.4°) 

And, for the five WMAP frequencies alone, there are: 

- minima at cosec \b\ = 2.15. 

In the North alone at the smallest angle (cosec \b\ = 4, 
ic \b\ = 14.5°) there is an excess over the 'best straight 
line' growing systematically with frequency, reaching 17% 
at 93 GHz. There is no equivalent in the South. 

Such fine structure, even when the data are integrated 
over all longitudes, lead us to believe that significant varia- 
tions, correlated in frequency, undoubtedly exist for smaller 
longitude bins - as we claim. It is very doubtful if such 
effects wo uld be subtracted out in the WMAP cle aning 
procedure ijBennett et alJl2003ol : iKomatsu et al]l2003fo - in 
which the weights vary considerably from one frequency 
to another. (They are, in order of increasing frequency, 
23-93 GHz, +0.109, -0.684, -0.096, +1.921 and -0.250). 
Furthermore, the weights are the same for N and S, whereas 
the dependence of slope of AT vs cosec |6| is different for 
N and S. Specifically, the slopes are, for 23 — 93 GHz: 
(<5AT/<5(cosecjZ>j)) s/N = 1.50, 1.47, 1.67, 1.61 and 1.50. That 
for 408 MHz is lower: 1.10 although, because the dependence 
of T on cosec \b\ is less linear, there is an appreciable error. 



2.5 The Situation at High Latitudes 

It is instructive to examine the situation beyond |6| = 45° 
in some detail because this region had been expected to be 
far enough from the Galactic Plan e to be completely free of 
foreground jBennett et alJl2003tJ). 

Other work jHangeriet^ilTl20o3) , mentioned already, 
has suggested that the axis of symmetry of a large scale 
modulation of the CMB signal is in fact, not along the Galac- 
tic N-S axis but nearer the axis throu gh the Ecliptic Poles. 
These results will be discussed later in IWibig fc Wolfendald 
2005); here, we examine this aspect from a different stand- 
point. Ideally, with a large foreground contribution, we 
might expect the 'N-S axis' to be defined by the lowest 
signals in the N and S hemispheres. However, even in the 
case of conventional Galactic radiation signatures (eg 21 cm 
line of the HI column density, synchrotron radiation, etc) 
the minima are away from the Galactic Poles. Furthermore, 
as we have noted, some of our CMB parameters increase 
with latitude above |6| — 45°. Thus, we determine the di- 
rections of both minima and maxima for each hemisphere. 
The parameters studied comprise those from the CMB anal- 
ysis already given, viz ATcmb, W(£=10) and the correlation 
coefficient for CMB temperature with respect to the CR in- 
tensity (to be described). In addition, we include directions 
for low (a nd high) spectral exp onents for CR protons and 
electrons jFathoohi et al.lll985 a b), making 5 parameters in 
all. With the exception of p and e, the directions are im- 
precise in that they are derived from the Quadrant - wide 
regions and |6| ranges 45° — 60° and > 60°; poor statistics 
preclude a more precise study. 

The results for the range of Z and b, and the mean values, 
are as follows. 

Maxima 

N: Z-range: 120°, mean 2:90° ±30° 
6-range: 20°, mean b: +50° ±5° 



S: Z-range: 90°, mean Z:230° ± 25° 
6-range: 15°, mean b: -50° ±5° 

Minima 

N: Z-range: 180°, mean I: 200 ° ± 50° 

6-range: 20°, mean b: +55° ± 5° 
S: Z-range: 130°, mean I: 20° ± 35° 

6-range: 20°, mean ?>: —55° ±5° 

Inspection of the values shows that indeed the means 
of both maximum and minimum values do not centre on 
b = ±90°. 

The analysis is such that a uniform (ie random) set of 
directions would have a range of Z - 360° and a mean with 
an error in t of ±90°, ie bigger than observed. 

The result is that the directions are consistent with con- 
formity as expected if foreground effects are significant at 
the high latitudes - to the extent of about 2-3 sigma in the 
North and somewhat greater in the South. The small errors 
for the fe-values are perhaps illusory because of the small 
overall range of b considered: 'Afe = ±11.25° (from the cen- 
tre of one bin to the centre of the other). 

Interestingly, the South, with its higher power, Fig- 
ure Q is^rwre sigjuficant. Comparison with the other result 
jHansen et alj f2004 ) is given later. 

At this stage, we can claim some support for the hy- 
pothesis that the foreground effects contribute to 'high' lat- 
itudes although to use the statistical weights just referred 
to is unwise since we run the risk of double counting, having 
already used the e and p result and some aspects of AT - 
and power - variations. 

2.6 Comparison with other studies of 
asymmetries in the cleaned maps 

Mention has already been made of the 'Gaussianity' differ- 
ences. Work on comparing t he sta tistic for each Galactic 
hemisphere jEricksen. et alJ (l2004aT) , b)has shown signifi- 
cant differences from expectation for the North but not for 
the South. The nearest we come to an effect of this sort is 
the variations from bin to bin in latitude of the power per sr 
at low I (I ~ 10); Figure^] In two of the Quadrants, 3 and 4, 
the variation is bigger in the North than in the South. The 
conclusion that there is 'an unexpected power asymmetry 
between the N and S hemispheres in the WMAP data' is 
certainly supported by our work. The authors of the Gaus- 
sianity work referred to above have gone on to use both the 
power spectrum and N-point correlation functions to argue 
that the axis of maximum asymmetry is close to the ecliptic 
axis. A somewhat similar result has co me from an analysis 
of 'the local curvature of WMAP data' jHansen et alJl2004l) 
although here anomalies are found for both hemispheres, 
at least for smoothing and the maximum asymmetry is be- 
tween the ecliptic axis and the Galactic Pole. Interestingly, 
it is demonstrated that for 5° — 10° smoothing the Southern 
hemisphere yields results too close to the Gaussian expected 
values. Inspection of the map for 'lake counts' ('lakes' are 
regions where the differential of temperature vs position is 
changing from negative to positive) shows maximum fluc- 
tuations which correlate reasonably with the 'hot regions' 
in Figure El In fac t, inspection of the map presented in 
lHansen et alJ i2004h shows that the centroids of the N and S 
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hemisphere regions having unusually high fluctuation values 
are somewhat away from the celestial poles. 

Specifically, they are at I, b:90° , +50° and 270°, -45°. 

These values are close to our estimates for the mean 
positions of the maxima 

- I, b: 90° ± 30°, +50° ± 5° and 

- l,b: 230° ± 25°, -50° ± 5°. 

The near agreement can be regarded as a confirmatory 
check on our analysis - and thus an agreement - to some 
extent - with our contentions. 



2.7 Power spectra for 'red' and 'blue' maps 

Before moving to the possible 'cosmic-ray-connection' it is 
necessary to refer to a feature that we have discovered which 
may have significanc e. We define points on the WMAP 
llTegmark et al.ll2003l) having positive values as 'red' and 
negative values as 'blue'. Power spectra have been derived 
for the whole sky for red maps and blue maps with the result 
shown in Figure Figure shows the ratio. Figures |SJ; 
and [SJi show the results for \b\ > f0°. It is evident that 
in both cases, \b\ < 10° and |fc| > 10°, there is a system- 
atic difference in the shape of the power spectrum between 
observation and expectation. 

Comparison with the 'limits' for the artificial universes 
shows two main features: 

i) a clear excess of red over blue at the highest ^-values; 

ii) a less significant deficit for £ < 10. 

The former can immediately be attributed to the contribu- 
tion from discrete sources. The significance of the important 
deficit at low I needs discussion. 

Insofar as the power levels at I values exceeding A£ ~ 5 
are largely independent, there are two values (at I = 5 and 
10) in this vicinity which are in the region of the lower limit. 
If this limit corresponds to the one standard deviation level 
then the chance probability is ~ 7%. The significance for 
cosmology of the difference at low t will be considered later. 

Of greater importance is the observation of the clear 
slow increase in ratio with increasing values of t. This in- 
crease is due to increased 'foreground' in terms of both 
Galactic and Extragalactic ('discrete' sources) components. 
Inspection of artificial Universe results shows that the 
chance probability is less than about 3%. Taking low I- values 
as a datum, and assuming that the total foreground affects 
only the red points, the result is that for the important one- 
degree peak the extra contribution is about 10%. Thus, from 
this argument, the power of the peak should be reduced by 
this amount. 



3 COSMIC RAY - CMB CORRELATIONS 
3.1 General Remarks 

As remarked earlier, the presence of a cosec \b\ variation 
in the mean CMB temperature at e ach and every WMAP 
temperature from 22.8 to 93.6 GHz, iBennett et"al"l (12003d) 
shows that there is a definite contribution to the uncleaned 
maps from Galactic - 'thin disk'- mechanisms. These are 



generally regarded as: electron synchrotron radiation, free- 
free radiation and dust. Insofar as the first mentioned is due 
to cosmic ray electrons and the dust is heated, in part at 
least, by CR protons, we expected 'cosmic rays' in general 
to be implicated in the foreground problem. It is true that 
subtracting off a cosec \b\ dependence, and arranging that 
the mean CMB temperature is 2.74K, removes much of the 
CR-related signal but some inevitably remains. Presumably, 
the residuals remain on a variety of angular scales and a vari- 
ety of magnitudes. Insofar as our knowledge of the distribu- 
tion of cosmic ray intensity and gas density (of all densities, 
including that in dust clouds) - in 3-dimensions- is very lim- 
ited, a determination of the likely CR-related contribution 
is very difficult. Some general points can be made, however, 
as follows: 

i) Lines of sight where the CR-intensity (of relevant en- 
ergy) is low might be expected to have low foreground 
(CMB) (denoted (CMB)/). 

ii) Lines of sight where the column density of gas is low 
('Galactic chimneys' and other regions) and thus low CR 
interaction should have low (CMB) / . 

iii) Lines of sight crossing regions where the CR in- 
tensity is high (SNR shel l s) dErlvkin fc Wolfendalel 1200 ll : 
IWolfendale fc Zhand Il994l : lOsborne et alJ Il995l) would be 
expected to have high (CMB) f . 

Each of these regions will be considered in turn, as will 
an all-sky corre lation using the low e nergy gamma ray data 
from EGRET llTegmark et alJl2003h and a more energetic 
component for |fc| < 10°. 

3.2 Correlations in Steep CR-Spectra Regions 

From the maps of electron and proton (Fat hoohi et alJ 
Il985al lbl) spectral indices we have chosen regions covering 
about 10% of the sky area. The results for a smoothing of 15° 
FWHM were taken for both the proton and electron compo- 
nents in that we expect CR 'variations' to be on this order of 
scale, corresponding to typical SNR-induced large scale tur- 
bulence. We expect the contours for protons and electrons 
to be similar, but not quite the same (differences of energies 
for the initiating particles, electrons and protons, will cause 
differences in diffusion coefficients). An estimate of the dif- 
ference can be made, as follows. The mean electron energy 
is about 1 GeV for the electrons responsible for the detected 
gamma rays and about 10 GeV for the protons so that the 
proton and electron escaping from t heir accelerating SNR 
will b e separated by of order 100 pc jErlvkin fc Wolfendalel 
2001) after a time of 10 6 y (ten times the lifetime of the 
SNR). At a typical distance of 1 kpc the angular separa- 
tion will be of order 6°. The somewhat disparate patterns 
of Figure|K|are thus understandable. 

In order to choose the most reliable regions of steep 
spectra - and because we are not sure which of electrons 
and protons are responsible - we have identified regions of 
common steep spectra for both proton s and e lectrons in Fig- 
ure (see also lErlvkin fc Wolfendalel (120011) 1: these should 
correspond to the most likely regions of maximum effect for 
the CMB. By 'maximum effect', here, we mean lowest CMB 
temperature because steep spectra indicate low particle in- 
tensities at the higher energies where we expect the extra 
pseudo-CMB contribution to arise. 
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Table indicates the regions so identified and the cor- 
responding values of (CMB) It will be noted that all are 
negative. The significance of the deficits will be examined 
shortly. 

3.3 Correlations for Galactic Chimneys 

Returning to the presence of chimneys of low column density 
of HI, in which we have claimed steep spectra (and where 
we used limited latitude regions above) , an analysis has been 
made of the mean CMB temperature along four major di- 
rections - lines of constant longitude (of width ±6°) over 
their whole latitude range; these directions do not always 
coincide with those of steep spectra. 

The mean temperature has been determined for 10° < 
|6| < 60°, this being the nominal region over which the chim- 
neys can be identified. The results are also given in Tabled 

It can be seen that our expected reduction in mean 
CMB temperature is observed in 6 or 7 out of the 8 regions. 

These results are potenti ally important and they wi ll 
be examined in more detail. In I Wolfendale fc Zhansl <1994l 
where evidence was given for correlations of steep spectra 
for low energy gamma rays with the Galactic chimneys and 
low HI regions, a comparison was made between b+ (+30°) 
and b- (—30°). It was found that the b+ results were better, 
those for negative latitudes being inferior, particularly for 
/ = 57° and 130°. These results related to electrons. With 
higher energy gamma rays, however, where protons are the 
progenitors, at negative latitudes, I = 57° was poor, but 
130° was good. The conclusion is that the anomalous posi- 
tive result for b- at / = 57° can, perhaps, be understood, not 
least because inspection of Galaxy-wide HI column density 
maps shows that the 'chimney' at negative b and I ~ 57° is 
barely visible. 

3.4 Correlations for Other Low Hi-Column 
Density Regions 

Another approach, resembling the previous one, is to look 
at regions where the HI column density is low, irrespective 
of the CR intensity and away from the Hi-chimney. The 
reason is that it may be that it is the HI density itself that 
is important rather than the CR intensity as such. Again, 
the results are shown in Tabled 

3.5 Correlations in High CR Intensity Regions 

If, as seems possible, it is CR interacting with gas that 
gives rise to a pseudo-CMB contribution, then there should 
be an excess of AT (CMB) in regions where the CR 
intensity is high. Probably the best evidence for high 
CR intensities comes from studies of gamma rays (typ- 
ically above 100 MeV) from SNR shocks fehat et all 
1 19851: I Wolfendale fc Zhand Il993) . It has been known for 
some years that the famous Loop I SNR has somewhat 
highe r gamma ray emission tha n expected fehat et alJ 
Il985l: IWolfendale fc Zhansl Il994f) . Most came from the 
well-developed ridge which is seen in radio, eg at 408 
MHz, rising along b ~30° from the Galactic Plane. The 
gamma ray work just referred to also shows evidence for 
CR acceleration in the shell of another SNR, Loop III, 



jLandecker fc Wielebinskilll970l) albeit at a lower level. Yet 
another SRN, Loop II, seems to be present, but rather weak 
and is disregarded. 

Table shows the results for the leading edge of the 
Loop I SNR (The North Polar spur) , which is best identified 
in radio. Also shown are values pertaining to the regions at 
\b\ < 10° where the two Loops cross the Galactic Plane (Fig- 
ure [7| shows the relevant CMB temperatures). This region 
is chosen so as to maximize any potential signal (positive or 
negative) in view of the higher gas density here. 

It is reassuring to note that the mean values of AT are 
almost always positive. 

It remains in this section to consider the statistical sig- 
nificance of the correlation results so far. 

Figure Q shows a histogram for the AT values all over 
the sky. The abscissa is in standard deviations and relates 
to areas of the magnitude used in the analysis. The symbols 
above the plot are for the individual areas referred to above. 
It is evident that the divergences are, taken at their face 
value, very significant. Formally, the chance of getting such 
a situation is ~ 10~ 4 for the negative excursions and ~ 10 -3 
for those in the positive direction. 

It is, of course, possible to question the significance on 
a number of levels: why choose the particular bin widths of 
latitude, and the particular size of the low N (HI) regions? 
The answer is that they were chosen before the results were 
known. Nevertheless, it is best to conclude at this stage that 
the correlations are 'very suggestive'. 

3.6 Whole-sky correlations of gamma rays and 
CMB 

The EGRET data feunter et all Il997l : iTeemark et alJ 
2003)have been used to search for correlations over the whole 
Galaxy, although it is appreciated that the statistical accu- 
racy of the gamma ray data is not great. Because there is 
a strong latitude dependence of the gamma ray intensity 
we study the CR-CMB correlation in narrow (but increas- 
ing with l&l) latitude bins. We concentrate on the highest 
energy (1 GeV) for which data of reasonable accuracy are 
available for the whole Galaxy. 

The analysis has been made in two ways. Firstly, by tak- 
ing the data in large, Quadrant size, bins in longitude and 
successive latitude bins and by taking 10 degree radius re- 
gions round successive pixels in the sky. Both methods have 
their merits, and demerits. Starting the with former, Fig- 
ure|S|shows a map of the correlations in a manner such that 
the latitude - and Quadrant - dependence can be sought. It 
is evident that there are both. The Inner Galaxy has no over- 
all positive correlation, formally it is negative, but the Outer 
Galaxy has a positive correlation, particularly in Quadrant 
2. There is also a pronounced dependence on Galactic lat- 
itude. This behaviour is surely suggestive of a relationship 
between the gamma ray intensity and the average CMB tem- 
perature, although the interdependence of bins is a worry. 
It might be thought that one could take wider intervals of 
latitude, near the Galactic Places, but here the CR intensity 
is such a strong function of latitude that artifacts can easily 
occur. 

Before continuing the (largely negative) correlations in 
the inner Galaxy will be considered further. It is possi- 
ble, in principle, that the CMB radiation is absorbed in 
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a particular geometrical way by high temperature ionized 
gas in the Galactic Halo, but the density and temperature 
needed are excessive. Another source of CR-intensity might 
be the answer: Inverse Compton interactions of electrons on 
starlight. Another possibility might be associated with the 
fact fSect. EOjl that there is the reverse gradient of the CMB 
power at latitudes above 10°. 

We found Quadrant-dependent effects in Sect. 12.21 and 
it is necessary to study Figure|H|in some detail. Although the 
latitude regions are of unequal width they can be considered 
with equal weight because, due to the strong ^-dependence 
of the gamma ray intensity, there are similar numbers of 
detected gamma rays in each band. In terms of strength of 
the correlation coefficient, the net coefficients are: 

Quadrant I = 4.8; Quadrant 2 +60.4; Quadrant 3 
+27.3; Quadrant 4 -16.1; 

It is appreciated that adjacent latitude bins are not in- 
dependent but this dependence is not too large, in view of 
the observed quite wide spread in coefficients from one bin 
to another. Thus, we conclude that there is some evidence 
for a correlation in the Outer Galaxy. 

For the range above \b\ = 20°, the summed coefficients 

are: 

Quadrant 1 +4.7; Quadrant 2 —4.0; Quadrant 3 +4.3; 
Quadrant 4 —4.7. 

The differences here are evidently of no significance. 

The positive correlation, (as evinced by this technique - 
but see later), is confined to \b\ < 20° and is predominantly 
in the Outer Galaxy. To what extent this result militates 
against a CR-CMB correlation at \b\ > 20° is not clear, in 
view of the poor gamma ray statistics at these latitudes and 
the fact, mentioned already, that rather than the average 
CMB temperature that is important it is the temperature 
of the 'background' over large angular scales (ie omitting the 
hot 'blobs') that should be important. This topic is taken 
up again later. 

Returning to the whole sky, the North- South difference 
is also of interest. The net coefficient from the North (pos- 
itive b) is +28.0 and that from the South is +38.7; again 
an S-excess. Unfortunately, it is not clear how to work out 
the statistical significance of this result, except to say that 
it is in the correct sense. What can be done, however, is to 
rotate the CMB map and rework the correlations. 

This has been done and we have gone further to the 
extent of examining the correlation between the CR-CMB 
correlation coefficient and AT, this has been done because 
we expect the coefficient to be higher in the regions of posi- 
tive excess temperature if, as is expected, CR (-like) effects 
give additional signals. 

The latitude region can be divided into two ranges, as 
usual: \b\ < 10° and |6| > 10°. For the former, all 3-latitude 
bands show significant positive signals (averaged over AT), 
for both N and S. Furthermore, the profile of the correlation 
of the coefficient, C(AT), with AT is similar in the N and 
the S. In the rotated maps the similarity between N and 
S is usually only strong in the region \b\ < 2°, and occurs 
because of the correlation between adjacent bands - it will 
clearly be a maximum for such a narrow band, of width only 
2°. Taking the average amplitude alone, the chance of seeing 
positive signals of the magnitude seen or greater is estimated 
to be < 1%. The shape similarity (or coefficient vs AT in the 



N and the S) will reduce this probability somewhat further, 
as will be demonstrated. 

Continuing to higher latitudes, \b\ — 20° — 30° and 
\b\ > 60° show small positive average coefficients whereas 
j&j = 30° to 45° and \b\ = 45° - 60° show small negative 
correlations. Before writing off evidence for a genuine CR- 
CMB correltion coefficient vs temperature at |6| > 20°, from 
this particular analysis, however, we draw attention to three 
features: 

i) In the S region (our preferred hemisphere for extra con- 
tributions) there is a distinct increase in correlation coeffi- 
cient with increasing AT in each of the 4 \b\ ranges. In the 
rotated maps this happens in only 20% of the individual 
cases and in less than 10% for each in a set of 4. 

ii) From past arguments, we expect there to be a corre- 
lation between the CR intensity and the 'power per unit 
solid angle' at small I. The average temperature has a bigger 
contribution from the high temperature, spatially condensed 
regions, regions which we do not expect to be affected - ex- 
cept in their base level - by CR (-like) effects. Taking all 
data at \b\ > 20° we have examined the correlation of the 
CR-CMB correlation coefficient with the mean power for the 
bin (Figure|3J 1=10). The coefficient is positive at the 15% 
level. 

iii) The behaviour above |fo|=45° has relevance to the situ- 
ation f or ATcmb referred to e arlier and to the work consid- 
ered in IWolfendale fc Zhansl Jl994l) . viz the presence of an 
axis of symmetry; inspection of the correlation coefficients 
(Figure|HJ for \b\ = 45° shows pseudo- Galactic N and S poles 
at: l,b= 80°, +55°; 220°, -55°. These directions are those 
of the maxima in the two hemispheres beyond \b\ = 45°. 
They are reassuringly close to those for AT. The uncertainty 
in angle is, again, approximately ±25°. 

The second approach, referred to above, has been made, 
using smaller spatial regions for the correlations. This in- 
volves comparing the appearance of the genuine WMAP and 
those for artifical universes. Two features in Figure |S] arc 
considered by us to provide the evidence for a bona fide cor- 
relation of CR with CMB: the excess in the Galactic Plane 
and the presence of symmetry about the plane. 

For the genuine map we derive the integrated length of 
longitude for which the CMB temperatures are above the 
'red' minimum (AT > 300,u.ft:). This is n 5° for the WMAP 
and, for the artificial universe, the mean is 33° ± 17°. Thus, 
formally, the WMAP effect is at 5 standard deviations. The 
significance is almost certainly less than that for a Gaussian 
distribution, but is it probably greater than corresponding 
to a 1% probability of being spurious. For symmetry we 
count the number of symmetrical regions of length at least 
Al — 10°. There are 7 symmetrical regions compared with 
2.0 ± 1.4 for the artificial universes,ie a 3.6 standard devia- 
tion effect. The corresponding chance probability is almost 
certainly less than 5%. 

Continuing to higher latitudes, we find, at |&|=30° and 
60°, integrated red regions significant to the extent of 2.8 
sigma and 2.6 sigma in comparison with the artificial uni- 
verse results. 

Our conclusion from this large scale CR-CMB correla- 
tion vs 'temperature' study is that it exists. 



10 T. Wibig and A. W. Wolfendale 



3.7 Restricted region studies at the highest 
energies 

EGRET data are available to 30 GeV for a restricted region, 
|6| < 10°, and a correlation analysis has been made here. 
Again, because of the strong 6-dependence of the gamma 
ray intensity, we take narrow bands of latitude. The results 
are shown in Figure [§] Although the spread is large there is 
a clear correlation of the high energy EGRET flux and the 
CMB excursions, for the region away from the very lowest 
Galactic latitudes, viz \b\ = 60° — 10°. 

We would not have expected correlations for |fc| < 2° 
(and to a lesser extent for 2° — 6°) where the lines of sight 
are very long and where there is considerable 'confusion'. 
For 6° < |6| < 10°, where we expect the best correlation 
following the arguments contributing to Figure |7] the slope 
is 2.2 standard deviations from zero taking all the data and 
significantly more for the (CR-favoured) S hemisphere, with 
smaller significance for the North. 

3.8 Indirect studies of a CR-CMB Correlation for 

\b\ > 10° 

Circumstantial evidence comes from N-S, fe-dependence ef- 
fect and Quadrant differences. 

For CR-related effects we would expect: 

A southern excess; mean for Q2 > mean for Q3; a fall-off 
with increasing latitude, where the quantity under examina- 
tion is either the power at low £ or the mean temperature. 

The results are, for Ql, 2 and 3 (Q4 is omitted, because 
of the 'big blue hole') 

1. Power at £=10. S/N>1, Q2<Q3, Figure fall-off with 
increasing \b\ (figure 

2. Mean CMB temperature. Southern excess S-N = 30[iK, 
similar for both Q2 and Q3. Q2 - Q3 = +60fj,K. A fall-off 
of (T) with increasing |6| from as far as \b\ = 45° (from right 
down to \b\ = 2° - 6°). 

Evidently, most of the features expected do occur, the 
only exceptions being the reversal of Q2>Q3 seen for the 
£=10 power intensity and the termination of the fall in (T) 
with increasing \b\ at 45°. The former is considered later in 
Sect. [I] The last mentioned is interesting in its own right. 
It is seen in all Quadrants (Figure and usually for both 
N and S. It is germane to remark that the 408 MHz and 5 
individual WMAP frequency plots all show an upturn (or 
constancy) of intensity at |6| = 60°, in order of increasing 
|6|. In our case, for the N, the increase (for Q's 1,2 and 3) is 
51 [iK whereas for the South it is +3/iisT; in the 408 MHz 
and WMAP frequencies, 4 out of the 5 sets of data which 
show an increase have a bigger effect in the North. In the 
power spectrum, £=10, data (Figure there is an increase 
in 4 out of 6 cases although that for the South exceeds that 
from the North. 

3.9 Summary of the CR-related - CMB 
correlations 

It is concluded that for |6| < 10° there is good evidence for 
a correlation of some form; each and every analysis shows 
an effect. The overall chance probability is less than about 
2%. 



In the range |6| = 10° to 20° there is ambiguity, epit- 
omized by the systematically low intensities in the power 
spectrum for both £ = 10 and £=100. This is reflected in 
(AT)only if Q4 is included. Without it, the mean tempera- 
ture for |6| = 10° — 20° is only a little above zero. 

Moving to |6| > 20° there is evidence for a correlation 
of CR with either power intensity (at £ = 10) or (AT) but 
the most direct evidence comes from the deficit in (AT) in 
the special regions (Galactic chimneys and regions of steep 
CR spectrum) and the excess in the SNR shock regions. 
Insofar as the special regions cover the range \b\ — 10° — 20°, 
too, there is evidence for correlations here, too, despite the 
ambiguities referred to above. The overall chance probability 
is, conservatively, at the 2% level. 



4 COMBINATION OF EXTERNAL AND 
INTERNAL STUDIES OF THE WMAP 

By 'external' we mean the comparison of the WMAP data 
with cosmic ray parameters and by 'internal' we mean a 
comparison of WMAP data within itself, eg N-S asymmetry, 
^-dependence. 

The largest effect concerning non-uniformity relates to 
low latitudes, typically \b\ < 10°, or somewhat wider. Both 
internally and externally there are differences at low lati- 
tudes; a Galactic origin is surely very likely here which may, 
or may not, be CR-related. The main question to be ad- 
dressed is "are there foreground effects beyond |6| = 10° or, 
in view of the anomalous region \b\ = 10° — 20°, beyond 20° 

This aspect will now be addressed in some detail. 

The evidence favouring CR-related effects at \b\ > 
20°can be summarised as follows: 

i) The deficit in the mean WMAP temperature associated 
with the regions of steep CR spectra and low HI column 
densities are quite significant, and extend to |6| fa 60°. 

ii) The excess for the North Polar Spur similarly extends 
to b = +60°. It, too, is significant. 

iii) The large scale CR-CMB correlation adds a little weight 
- at least at the 3% significance level (although the absence 
of much correlation in the Inner Galaxy is a worry). 

iv) The N-S symmetry is large at small I (for the power 
spectrum) and continues to \b\ ~ 55° — 60°, after which 
it falls. This is mirrored by the 408 MHz signal and by the 
individual WMAP data for the various frequencies. It is gen- 
erally agreed that the |6| -dependence of the 408 MHz and 
the individual WMAP data sets is due to CR-related Galac- 
tic foreground and the low I power spectrum asymmetry 
presumably has the same cause. 

v) The form of the actual latitude dependence of various 
parameters presumably also has the same cause. For the five 
frequencies the ratio of AT(20° - 30°) to AT(90°) is 2.1 (it 
falls slowly with increasing frequency from 2.3 to 1.9). Our 
ratio for the power at I = 10 is 1.7, ie 1.3 for the root 
mean square temperature. The ratio for larger ^-values will 
be smaller. The relationship between the longitude-averaged 
AT values for the five frequencies and that for V AT 2 for 
small £ (I = 10) is not completely clear but if the AT fluc- 
tuations are a constant fraction of the total DC signal then 
the '2.1' is to be compared with the '1.3'. Taken at its face 
value it would mean that the maps were cleaned to the ex- 
tent of a factor 4 - and a greater factor at higher I (eg by ~ 8 
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at £ = 30, using Figure[TJi). These factors relate to |fc| > 20° 
where most of the CMB signal used for cosmological anal- 
yses resides. A lack of cleaning to the extent derived (25% 
of the signal at £ = 10 being CR-related noise) at I ~ 10 is 
clearly serious. 

vi) At the highest latitude the consistency of the maximum 
intensity, from one parameter to another, adds confidence. 

vii) A comparison of both (T) and power (for I ~ 10) for 
Q2 and Q3 gives understandable results in terms of Galactic 
phenomena. 



5 APPLICATION TO COSMOLOGY 

5.1 General Remarks 

In the absence of clear evidence as to the exact origin of the 
foreground responsible for the various anomalies in the CMB 
map it is not possible to correct for it and make an updated 
cosmological interpretation. However, some progress can be 
made by looking in detail at three important parameters: 

- the height of the 'one-degree peak'; 

- the lvalue of the peak, and 

- the amplitude of the power spectrum at low I (£ < 10). 

5.2 The height of the one-degree peak 

Figure tTUI shows the power spectrum for the whole data used 
by us. Also shown is what is generally regarded as the stan- 
dard form as given by the WMAP workers themselves (and 
a curve 10% below this). We note that 'our' value is low. 
Interestingly it is nearer to the conclusions drawn from pre- 
WMAP observations ( see iKomatsu et al.1 <2003h 1. Specifi- 
cally, the heights of t(£ + l)C t /2i:{^K 2 ) are: 

- WMAP group iBennett et alJl2003b[) 5700 

- previous work [ jKomatsu et alll2003h ] 4925 ± 175 

- present work 4250 (whole sky) 3490 (|6| > 10°, scaled 
up). 

The difference i n height of the p eak b etween our value 
using the data of [Teemark et alJ (|2003)and the datum 
IBennett etafl (2003b) is seen to be at least 25% and there- 
fore serious. 

Turning to the variations from Quadrant to Quadrant 
and N to S, Figures Q and indicate that these are (at 
I ~ 200) of order 10%; similarly, Figure [£] suggests that the 
observed peak is some 10% too high, in the sense that the 
red (with extra signals) and blue (without) differ such as to 
give this difference. 

There are two ways forward at this juncture, the first 
is to assume that our derived value of 3490 (less 10%), 
ie 3140 \xK 2 is correct and the other that there is a sys- 
tematic displacement i n the scale of the resu lt derived us- 
ing the cleaned data of iTeemark et all (I2003T) . and that we 
should t ake 90% of the WMAP groups' result IBennett et alJ 
j2003bMHinshaw et alJ d2003l) . viz 0.9 x 5700 = hUO^K 2 . 

The implication of a low peak value is most likely that 
Cl m , the mass content of the universe, should be reduced, 
although we appreciate that other variants to the cosmo- 
logical parameters can be taken. Nevertheless, a reduced 
peak height is mos t easily a chieved by increasing Q. m (or 
specifically n m h 2 ) llHull2003f) . A peak height of 3590 /j,K 2 



would indicate Q. m h 2 — 0.40, a very high value, and 5130 
fj,K 2 yields Q. m h 2 = 0.20, itself significantly higher than the 
'standard value'. 

5.3 The I- value of the one-degree peak 

Less susceptible to systematic errors is the lvalue of the 
peak. Conventionally it is at £a = 220 but it is apparent 
that our result is systematically lower. 

Values have been derived for \b\ > 10° for the Quad- 
rants, by fitting a 'standard shape' of W(£), derived from 
the whole data. The values are: 

- Quadrant 1 205 ± 10; 

- Quadrant 2 180 ± 10; 

- Quadrant 3 190 ± 10; 

- Quadrant 4 210 ± 10 

giving a best-estimate of 196 ± 5. This value is 24 ± 5 lower 
than the conv entional 2 20. Using the expression A£a/£a = 
-0.24 AQ m h 2 llHull2003!) and a datum of Q, m h 2 = 0.15, we 
have AQ, m h 2 = 0.17, ie fi m /i 2 =0.22. 

An alternative approach is by way of Figure |5] Again, 
using the blue spectrum we have £a = 200 ± 10, ie A£a = 
-20±10 and AQ m h 2 = 0.06, ie Q m h 2 = 0.21. The estimated 
error is ±0.05. 

5.4 The situation at low f-values 

It is well know that there is a loss of power below I = 10 
in com parison with conventional expectation iBennett et alJ 
l2003bl) and this has given rise to speculation of 'new physics' 
(eg Efstatiou (2003)). One possibility is that it is 'mere' cos- 
mic variance that is responsible and it is true that this is 
large at low £- values. The big differences from Quadrant to 
Quadrant, both for low and high latitudes, have relevance 
here; they show that there may be genuine differences. How- 
ever, these may be due to CR-related foregrounds insofar as 
the deviations from one Quadrant to another are bigger in 
the region |6| < 10°. 

At first sight there is evidence for an even greater fall 
off of power at low £ with increasing latitude so that, in 
the absence of a cosmic- variance-inspired effect, the anomaly 
might be thought to be secure. However, we are struck by 
the results in Figure |S] which shows a near-normal power 
law intensity at low £-ior the blue regions, ie the favoured 
set of data. Although the peak of the blue curve is low, from 
£ = 100 downwards, the form of the blue curve mirrors that 
of the conventional prediction all the way down to the lowest 
values of £. A possibility is that the form of the foreground 
contribution is such as to reduce the apparent power at low 
£-; we remember that the total temperature variations are 
normalized to zero. 



6 CONCLUSIONS 

The major conclusion from this work is that there are fea- 
tures of the WMAP of the fluctuations in the CMB that 
are not related to the early universe, but are due to a 'fore- 
ground'. Some of the foreground at least is surely of Galactic 
origin and it appears to exist at all Galactic latitudes. 

Large scale asymmetries (also seen by others) exist on 
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scales of North vs South and Quadrant to Quadrant. Our 
understanding of them is limited but guidance seems to come 
from correlations of the mean temperature, and intensity of 
the power spectrum, with cosmic ray features. 

The implications for cosmology are complex, but some 
remarks can be made, as follows: 

1. Presently derived cosmological constants are more un- 
certain than usually appreciated. 

2. If it is assumed that the only parameter to be adjusted 
is the universal mass fraction Q. m h 2 then there is evidence for 
it being larger than the canonical 0.15. Two of our analyses 
suggest that it is ~0.22 ±0.05. Another analysis gives a value 
as high as 0.40. It is premature, however, to put much weight 
on this very high value. Nevertheless, a value higher than the 
canonical Q. m h 2 — 0.15 is preferred. Of course, variations of 
the values of other cosmological parameters are possible, but 
such variations would also be of great interest. 

3. The situation at low £ is interesting and our analysis 
suggests that the more conventional view may prevail, viz 
that the apparent missing power may be an artifact. 
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Figure 1. Power spectrum of the WMAP data iTegmark et alJl2003Tl for the two Galactic hemispheres (a) North and (b) South. Moving 
upwards, the latitude ranges are (plus or minus) : < 20°, 20°/30°, 30°/45°, 45°/60° and > 60°. The data for positive latitudes (the 
North) have been smoothed and also plotted in (b) to allow observation of the fact that the 'Southern Excess' continues to 'high' 
latitudes. Successive spectra have been displaced upwards by a factor 10 to aid appreciation, (c) Power spectrum South/North ratio vs 
I for all the data; the 'cosmic variance' is also shown, (d) South/North ratio for the various latitude bands (there is a smoothing due to 
interdependence at low latitudes). 
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Figure 2. (a) Power spectrum, for latitudes below 10°, for each Galactic Quadrant in turn. The lowest line refers to Quadrant 1, and 
so on. Successive lines are displaced upwards by a decade to help appreciation of their shapes, (b) As (a) but for latitudes above 10° . (c) 
The ratio of the power in the other Quadrants to that in Quadrant 1, for |b| < 10°. (d) As (c) but for |b| > 1 0° . Overall Power spectrum 
and spectra for the two latitude ranges for (e) Quadrant 1 and (f) for Quadrant 3. Also shown is the result I Wibig & Wolfendale 2004) 
for the whole sky divided by 4 to give the expected power per quadrant. 
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Figure 3. Mean CMB temperature for each Quadrant for North (solid lines) and South (dashed lines) (a-d), and average values for the 
whole Quadrants (e,f) as a function of Glactic latitude. 
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Figure 4. (a) and (b) Power at i = 10 for the different latitudes and Quadrants, (c) and (d) as (a) and (b) but for i = 100. The power 
is expressed per unit solid angle. The open circles relate to the S-hemisphere and the filled circles relate to the North. It is evident that 
there are significant variations. 
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Figure 5. Power spectrum for red (positive) and blue (negative regions) and their ratios. Studies of artifical universes show essentially 
no red/blue differences. show the limits 
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Figure 6. Map of the sky (Galactic Centre at the centre, longitu de increasing to the l eft) showing regions of steep spectra for protons 
and electrons I Fathoohi ct al. 1985a) and minima in the WMAP jTegmark et alJl200.'t . The thick lines represent low column densities 
of atomic hydrogen (chimneys) and interarm regions, 'e' represents two regions of steep electron spectra below the low CMB region. The 
ridge ('North Polar spur') for the Loop I SNR is also shown as a dashed line. 
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Figure 7. Excesses and Deficits. Values from Table HI for the C-R associated excesses and deficits in comparison with the frequency 
distribution of mean temperatures for similar spatial regions shown by solid histogram (dashed and dotted histograms are for Noth and 
South hemispheres, respectively). 
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Figure 8. Gamma rays above 1 GeV f lHunter et aL^lSST!) and see iWibig fc Wolfenrtalfll2*'o04l) ) vs CMB temperature iTegmark et al.l 
2003). Map of the correlation coefficients with the separate Quadrants identified (Galactic scale as in Fig.|HJ. The units of the correlation 
coefficients are 10" 2 . The latitude ranges arc: |6| < 2° , 2° - 6° , 6° - 10°, 10° - 20°, 20° - 30°, 30° - 45°, 45° - 60° and above 60°, with 
similar divisions (below b = —2°) for negative latitudes. Although the statistical precision of any individual coefficient is poor there are 
some important patterns, most notably that there is a stronger correlation in the Outer Galaxy, particularly in Quadrant 2 (in which 
all astronomical indicators are stronger), and that there is a concentration in the Galactic Plane. We think that the large gamma ray 
contribution from Inverse Compton interactions and other factors, may have spoiled the correlation in the Inner Galaxy. 
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Figure 9. Gamma rays at 30 GeV from iHunter et al.lll997t) for the restricted latitude range where adequate statistic precision is 
available, viz |b| < 10° vs CMB mean temperature iTegmarlTet al.l20 031 The best-fit lines show that the correlation is rather strong for 
the range |6| = 6° — 10°. Open circles: Southern hemisphere, closed circles: Northern hemisphere. The large circles represent the data 
concentrated in 3 bins. 
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Figure 10. The upper curve is a standard result from WMAP iBennett et aljfeoOSat LSpereel et"afll200a|l . t he next one down is 10% 
below this and the next p air are from our own analysis of the WMAP data for the whole sky presented in iHunter et al.lfl997f) (see 
iWibig fc WolfendaldEool '). 



